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OO ' Abstract 

o . 

■ We address the question of which phase space functionals might represent a quantum state. 
I We derive necessary and sufficient conditions for both pure and mixed phase space quantum 

Tiij- ' states. From the pure state quantum condition we obtain a formula for the momentum cor- 

■ relations of arbitrary order and derive explicit expressions for the wavefunctions in terms of 
time dependent and independent Wigner functions. We show that the pure state quantum 
condition is preserved by the Moyal (but not by the classical Liouville) time evolution and is 
consistent with a generic stargenvalue equation. As a by-product Baker's converse construc- 

r~| , tion is generalized both to an arbitrary stargenvalue equation, associated to a generic phase 
space symbol, as well as to the time dependent case. These results are properly extended to the 
mixed state quantum condition, which is proved to imply the Heisenberg uncertainty relations. 
Globally, this formalism yields the complete characterization of the kinematical structure of 
I Wigner quantum mechanics. The previous results are then succinctly generalized for various 
quasi-distributions. Finally, the formalism is illustrated through the simple examples of the 
harmonic oscillator and the free Gaussian wave packet. As a by-product, we obtain in the 
former example an integral representation of the Hermite polynomials. 



X 



1 Introduction 

The Weyl- Wigner approach to quantum mechanics [1]-[16] constitutes an alternative consistent 
formulation of quantum theory. Its main advantage resides in the fact that both observables and 
states are represented by ordinary c-numbers in a non-commutative phase space and that this 
representation strongly mimics classical statistical mechanics. To any given quantum operator A 
from the quantum algebra A we can attribute a " classical" phase space counterpart A in the classical 
algebra A through Weyl's correspondence rule [2]. Likewise, there is a phase space entity called the 
Wigner function F^{x,p) [4], which plays the role of the density matrix. Weyl's correspondence 
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rule endows the phase space with the structures of a non-commutative product and a bracket (the 
Moyal bracket). All the usual quantities of interest in operator quantum mechanics (probabilities, 
mean values) can be evaluated in phase space through simple formulae reminiscent of classical 
statistical mechanics. 

The intuitive methods of deformation quantization have found interesting applications in a wide 
range of fields of research. Standard applications in non-relativistic quantum mechanics include 
the fields of the semiclassical limit of quantum mechanics [1, 9, 17, 18, 19], hybrid dynamics [20], 
collision theory [1, 13] and quantum chaos [1, 21, 22, 23]. The deformation techniques have also 
been used to provide an alternative phase space formulation of quantum field theory [24, 25] and to 
address the quantization of constrained dynamical systems [26] . Moreover, the deformed algebraic 
structures play a key part in some current developments in M-theory [27] -[30]. The vast number 
of applications reinforces the general importance of the Wigner theory not only as an alternative 
formulation of quantum mechanics but also as a comparatively simpler formalism where to test and 
further develop the deformation approach. 

A problem lying at the core of Wigner theory is whether its formulation is fully consistent with 
standard operator quantum mechanics. The logical equivalence of the two formulations presupposes 
that we should be able to formulate quantum mechanics in phase space autonomously and subse- 
quently be able to rcdcrive operator quantum mechanics. This program has had a long history. In 
refs.[2]-[9] the phase space formulation of quantum mechanics was firmly established. In refs.[8] and 
[10]-[12] the inverse route was explored: starting from phase space one tried to recover the operator 
formulation. In our opinion this path has not been fully completed. 

In this work we wish to reinforce the general belief that the phase space formulation of quantum 
mechanics is a self-consistent autonomous theory in its own right by addressing the following issues: 

(1) given some real, normalizablc phase space function F^{x,p) how can we assess whether it is 
a bona fide Wigner function?; (ii) how can we tell whether F^{x,p) represents a pure or a mixed 
state?; and finally, (iii) can we reconstruct the wavefunction or the density matrix from the Wigner 
function? Partial answers to some of these questions have been known for some time. We will 
review the known results and state what is novel in due course. 

The starting point of this paper is the presentation of several quantum conditions determining 
whether a given phase space functional represents a pure or mixed quantum state. These conditions 
provide a complete classification of phase space functionals in terms of: (i) non-quantum states, (ii) 
pure quantum states and (iii) mixed quantum states. One of the pure state quantum conditions 
presented was originally introduced by Tatarskii [12]. For completion we shall review his derivation. 
The rest of the paper is then dedicated to: (1) proving the compatibility between the quantum 
conditions and the main structures of Wigner quantum mechanics and to assess the degree up to 
which these structures are in fact determined by the quantum conditions. Namely, we shall prove 
that: (la) Moyal (but not classical Liouville) dynamics preserves the pure and mixed quantum 
conditions, (lb) the solution of a generic stargenvalue equation is a pure state and (Ic) the mixed 
state condition implies the Heisenberg uncertainty relations for general non-commuting observables. 

(2) Rederiving standard operator quantum mechanics from the phase space formulation. This 
amounts to: (2a) deriving explicit expressions for wave functions (or density matrixes) associated 
to a generic pure or mixed state Wigner function, (2b) proving that these wave functions satisfy the 
proper eigenvalue (Schrodinger) equation if the original Wigner function satisfies the corresponding 
stargenvalue (Moyal) equation. The existence of such associated wave functions was established 
[10, 11] for the energy stargenvalue equation. Here we fully generalize this result by considering 
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the case of arbitrary stargenvalue equations (associated to a generic phase space symbol) and the 
time dependent case. Moreover, and most importantly, we provide the complete specification of the 
associated wave functions and density matrix. Finally, (3) we extend the previous results to other 
phase space quasi-distributions stemming from alternative ordering prescriptions. 

This set of results provides the missing mathematical structure concerning the relation between 
standard operator and Wigner quantum mechanics and firmly establishes the equivalence of the two 
formulations. It also yields a (still missing) complete characterization of the kinematical structure 
of phase space quantum mechanics. This structure is not compatible with classical statistical 
mechanics. In particular the Liouville time evolution and the classical measurement collapse do 
not preserve the pure state quantum condition. On the other hand the Moyal theory (dynamics 
and stargenvalue equations) is fully compatible with the quantum conditions. Whether or not 
the quantum dynamics and the quantum measurement collapse are uniquely determined by the 
quantum kinematical conditions is an open question to be addressed in a future work [31]. An 
affirmative answer to this question will imply that quantum mechanics is the only logical statistical 
formulation for the constrained set of phase space distributions that satisfy the pure state quantum 
condition. 

This paper is organized as follows. In section 2 we briefly revise some notions concerning pure 
and mixed states in ordinary quantum mechanics, which we regard as pertinent for the sequel. In 
section 3 we present a summary of quantum mechanics in phase space. In section 4, we define 
pure states in phase space and explore some of their properties. In section 5, we derive and study 
the quantum condition for mixed states. In section 6, we briefly generalize the results to other 
quasi-distributions. Finally, in sections 7 and 8 we consider two examples: the time independent 
simple harmonic oscillator and the free Gaussian wave packet. As a by-product we obtain a new 
(to the best of our knowledge) integral representation for the Hermite polynomials. 

Finally, let us make the following remark: for simphcity we shall only consider the case of one- 
dimensional systems displaying a classical phase space with the structure of a flat manifold. The 
generalization to higher dimensional systems is straightforward but the extension to curved phase 
space dynamical systems might only be possible under some restrictive assumptions and, even in 
that case, may constitute a rather non-trivial problem [31]. 

2 Quantum states 

Let us then consider a one-dimensional system with fundamental observables g, p satisfying Heisen- 
berg's algebra: 

[q.pl^ih. (1) 

The state of the system is completely specified by a vector > in the Hilbert space or by the 
so-called density matrix p = \ip >< ip\. Let >, a; G 3?} be the set of eigenstates of q. The 
wavefunction is defined by iIj{x) =< x\ip >. If ip^x) is square integrable and normalized (with 
respect to the L^(K) norm) then it is an admissible pure quantum state. We interpret \ip{x)\'^dx as 
the probability of finding the particle in the interval ] x, x+dx [ . Likewise, \ip{p)\'^dp is the probability 
that the particle's momentum belong to the interval ]p,p + dp [. Here ■0(p) is the Fourier transform 
of ip{x). Planck's constant makes its appearance in eq.(l). The previous considerations then lead 
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to one of the main results of quantum mechanics, namely Heisenberg's uncertainty relations: 

Ax-Ap> ^, (2) 

where (A^;)^ =< > — < z 

Now, let A be some hermitian operator functional of q and p with discrete nondegenerate 
spectrum^. Suppose further that {|a >, a e /} is a complete orthonormal set of eigenvectors of A: 

A\a >= a\a >, a & I. (3) 

Then the wavefunction can be expanded as follows: 

|V'>=E^a|a>, (4) 

where {Ca, a e /} are a set of complex constants, such that Y^aei \Ca\'^ — 1- The probability that 
a measurement of the observable A yield the eigenvalue a is |Cap. 

Alternatively, suppose that instead of a single state >, we have a statistical mixture of states 
each one with probability pp. In this case the density matrix becomes: 

f3eB peB a,bel a,bel 

0<pp<l, V/3eS; = (6) 

Pes 



where 



and 



. Pa6 = E/36BP/3Cf)Cf*, a,bel. 

Prom these relations we can easily infer the following properties of p: 



(7) 



1) p is hermitian; 

2) p is non negative: 

< > > 0, (8) 

for any state |0 > in the Hilbert space; and 

3) p is normalized: 

Trp = 1. (9) 

Conversely, it can be shown [32, 33, 34] that if a matrix satisfies properties 1) - 3), then it admits 
the following decomposition: 

P = Pf^li^^''^ X ^^^^ (10) 

where ||-?/''^'^) >,P ^ i?| is some set of normalized states and pp are some real constants satisfying 
(6). Equation (10) means that the density matrix is inevitably a superposition of pure states. 

^The generalization to continuous and degenerate spectra is straightforward. 



4 



Conditions 1) - 3) allow for a probabilistic interpretation of the diagonal entries of the density 
matrix, which is independent of the complete orthonormal basis of the Hilbert space. 

It should also be noticed that we can always find a decomposition (10) where the states \'ip^'^'' > 
are orthogonal, i.e. < ijj^"'^ [i/j^^^ >= 6a,f3- In this case the density matrix is diagonal and the states 
> are eigenstates of p. The associated eigenvalues are non-negative and normahzed by 
virtue of (8) and (9). Using this particular decomposition one can easily prove that: 

i) A hermitian and normalized operator p (i.e. satisfying properties 1) and 3)) is a pure state 
density matrix iff it satisfies p = pp. In fact property 2) follows immediately from the preceding 
properties. Hence, p satisfies 1) to 3) and can be written in the form (10) where pp satisfies (6) and 
< >= Sa,i3- We then have: 

l3eB Pes 

Pp = for some /3o e -B, (11) 

where we used eq.(6). Hence, p — >< and so it is a pure state. 

We leave as a simple exercise to further prove that: 

ii) If p is a density matrix (i.e. it satisfies properties 1) to 3)) then p is a pure state iff Trpp — 1. 

iii) If p is a density matrix then it is a strictly mixed state iff Trpp < 1. 

This completes our review of some well-known results of the density matrix formulation of 
quantum mechanics. This section can then be summarized as follows. If p is a matrix which satisfies 
properties 1) - 3) then it is an admissible quantum state. If it further satisfies ii) (respectively iii)) 
then is a pure (strictly mixed) quantum state. Alternatively, if it satisfies 1),3) and i) then it is a 
pure quantum state. 



3 Weyl-Wigner formulation of quantum mechanics 

In this section we revise the Weyl-Wigner formulation of quantum mechanics. Let us start with 
Weyl's correspondence rule. The Weyl map V\v is an isomorphism from the quantum algebra A of 
linear operators with an operator product ■ and a commutator [, ] to the classical algebra A of phase 
space functions with a non-commutative ^-product (the Groenewold product [5]) and a bracket (the 
Moyal bracket [6]). It is defined by: 

A {q,p) ^ Vw{A) = A{x,p) ^ ^ j j dri Tr{A {q,p) e'€«+'''*^} g-'S^-^^f . (12) 

The definitions of the non-commutative *-product and the Moyal bracket follow immediately: 

A*B =Vw{A-B) 

[ABU ^^,{A*B-B*A) = ^Vw{ 
From equation (12), we get the following formulae: 

A*B = A(x,p)ef^B(x,p) = A{x,p - q£)B(x,p + i||) 

= ^/(ip'd/da;'da;M(a;^p')^(^^/)e^^^^'^'"''''^+^'^"''""^+^''^''"'^'^^ (14) 
JA,BU =^A{x,p)sm{^j)B{x,p), 
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where j is the " Poissori'' operator: 

(15) 

dx dp dp dx ' 

the derivatives d and d acting on A and B, respectively. The Weyl map of the density matrix is 
proportional to the Wigner function [4] : 



(16) 



In particular for a pure state, p = |^ >< we get: 

F^rei^,P, t)-^Jdy e-'y^r - ^, ^ + ^, ^ . (17) 

The pure state Wigner function is square integrable, normalized (in the sense that / dxdyF^^^{x, p, t) — 
1) and real, but, in general, it is not positive defined. The mean value of the operator A can be 
computed according to the following appealing formula: 

< A{t) j dx j dp A{x,p)F^{x,p,t). (18) 

If the operator A has nondegenerate spectrum and |a > is one of its eigenvectors (3), then the 
corresponding Wigner function F^{x,p) = ^%Vw >< a|) is the solution of the following *- 
genvalue equations [11, 16]: 

A{x,p) * F^{x,p) = aF^{x,p), F^{x,p) * A{x,p) = aF^{x,p). (19) 

The probability that a measurement of the observable A yield the eigenvalue a is then given by 

V{A{t) = a)= j dx j dp F^{x,p,t)F^{x,p), (20) 

provided F^ is normalized. If the spectrum is degenerate, then the generalization is straightforward, 
[16]. 

If the state is mixed, then from (5), we have: 

F'^(x,p,t) = Y.PpFZreA^^P^t) = E Pa,{t)FY,{x,p), (21) 
PeB a,bel 

where the non-diagonal Wigner function F^ reads: 

the pure state Wigner functions F^j.^j^{x,p,t) are of the form (17) (and can always be chosen so 
that they satisfy F^^^j^{x,p,t) * F^^^^^{x,p,t) = ^^F^^^j^{x,p,t)), pp satisfies (6) and pah is given 
by (7) thus satisfying the following conditions: i) pab = pla, h) Paa > 0, Va G / and iii) J2a&i Paa = 1- 
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If the wavefunctions ■^i, ■02 are solutions of the Schrodinger equation, 



then the non-diagonal Wigner function obeys the Moyal equation: 

dt 

with 



M' 



(24) 



H{x,p) = ^ + V{x). (25) 
Likewise if ipj are energy eigenstates with eigenvalues Ej {j = 1, 2), 

-|^^ + n^)^. = ^A, i = l,2, (26) 
then is a solution of the *-genvalue equations: 

H{x, p) * F^ = E^F^, F^ * H{x, p) = E^F^ . (27) 
This completes our overview of the Weyl- Wigner formulation of quantum mechanics. 

4 Pure states in quantum phase space 

In this section we will start to address the equivalence of the standard operator and Wigner's 
formulations of quantum mechanics. We shall first concentrate on pure states. Let us start by 
presenting the relevant definitions. A phase space function is a pure quantum state iff is of the form 
(17) for some complex valued, square integrable and normalized function ip{x,t). In other words, 
iff is the Weyl transform of a pure state density matrix. A phase space function of the form (22) for 
some complex valued, normalized functions ipai^) '4'b{x) will be designated by non- diagonal pure 
state. They will be important for the study of mixed states in the next section. When there is no risk 
of confusion we will also designate them by just pure states. Finally, we will also consider phase 
space functional of the form (22) where the associated wave functions are distributions. These 
states are non-physical (even if real the associated wave functions do not belong to the Hilbert 
space) and will be important (as in standard operator quantum mechanics) because they provide 
basis for the space of physical phase space functionals. We will designate them by distributional 
(non- diagonal, if complex) pure states. 

Let then F^ (x.p.t) be some function in quantum phase space. How do we know whether this 
function represents a pure quantum state? i.e. how do we know that this function is of the form 
(17)? The two following Lemmas provide two alternative characterizations of phase space pure 
states. 

Lemma 1: A real and normalized phase space function is a pure state Wigner function iff it satisfies: 



pw ^pw ^ _]_pw /28) 

pure pure 'ZttTi P"^^ V / 
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Proof: F^^g is a pure state Wigner function iff it is the Weyl transform of a pure state density 
matrix {ip >< ipl. By applying the Wcyl map to the conditions 1), 3) and i) of section 2 we 
immediately get that F^^^ is: 1) real, 2) normalized and 3) satisfies (28). 

Conversely, if a phase space function F satisfies (28) then by virtue of Theorem 2 (see ahead) 
F is of the form (22). If in addition F is real and normalized then F is of the form (17) for some 
normalized wave function ip. Hence, F is a pure state. □ 

Equation (28) provides an elegant, compact form for the pure state quantum condition. However, 
the evaluation of the starproduct in (28) may constitute in practice a rather difficult problem. 
Following Tatarskii [12] we now provide an alternative, more practical pure state quantum condition: 

Lemma 2: Let F^{x,p,t) be square integrable and let us define the following function: 

Z{x,j,t) = J dp e'^PF^{x,p,t). (29) 
If Z{x, j, t) satisfies the nonlinear partial differential equation, 

— lnZ(x,j,t)= y —\nZ{x,j,t), (30) 

then F^(x,p, t) is phase space function of the form: 

t)-^ldy ^''''^2 [x - ^, ?j i,, [x + y, , (31) 

where V'l, V'2 are some complex square integrable functions. If in addition F^^ is real and normalized 
then it is a pure state Wigner function. Conversely, if F^(a;,p, t) is a pure state (possibly non- 
diagonal) then it satisfies eq.(30). 



Proof: Eq.(30) can be seen as a wave equation, where the "velocity" is c = |, if j plays the role of 



h 

2' 

"time". Its solution is: 



In Z{x, i, t) = In ^2 (^x - ^, t j + In ^1 (^x + ^, t j , (32) 

where ipi, 1/J2 are arbitrary complex functions. On the other hand, Z(x,j, t) is the Fourier transform 
of F^{x,p,t). And so we get: 

F^(x,p, t)^^ldy e-'yPZ{x, y,t) ^ ^ J dy e^^^'V^^ - ^, ?j (^x + y, , (33) 

and so F^ is indeed of the form (31). Moreover and since F^ is square integrable so are ipi and 
ip2- If we further impose that besides being a solution of (30), F^{x,p,t) should also be real and 
normalized, then it is trivial to verify that F^{x,p,t) is indeed of the form (17), where ip{x) is 
normalized. Finally, the converse result follows immediately from the definition of Z[x,j,t).a 

We shall call eq.(30) the pure state quantum condition. Its real and normalized solutions are 
pure state Wigner functions. A solution of eq.(30) is a physical state iff is a pure state. The non- 
real solutions, on the other hand, are the non-diagonal pure states. Equation (30) was originally 
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introduced by Tatarskii [12]. Let us then explore some properties of the quantum pure states. 



Theorem 1: Let {x,p,t) satisfy the quantum condition (30) at the initial time t = and 
let its time evolution be given by the Moyal equation (24). Then F^{x,p,t) satisfies (30) for all 
times. 

Proof: This theorem was also originally stated by Tatarskii [12]. Let us briefly review his proof: 
if [x,p,0) satisfies eq.(30) then it is given by cq.(31) for some complex functions ipi{x,0) and 
ip2{x, 0). It is well known [6] that if ipi{x, t), ip2{x, t) are the Schrodinger time evolution of ipi{x, 0) 
and -02(2^,0) then the associated Wigner function F^{x,p,t) (given by eq.(31)) satisfies the Moyal 
equation for the Weyl symbol of the original Hamiltonian operator. Since the Moyal time evolution 
is uniquely determined by the initial time Wigner function F^{x,p, 0), we conclude that F^{x,p, t) 
is of the form (31) for all times. □ 

Later in this section we will show that this property is not shared by the classical (Liouville) 
time evolution. The next theorem extends Baker's converse construction to the solution of a generic 
stargenvalue equation. 

Theorem 2: Let A{x,p) = Vw{A) be a generic phase space symbol associated to a generic lin- 
ear operator A. Then the solution of the stargenvalue equations: 

A{x,p) * F^{x,p) = aF^{x,p), F^{x,p) * A{x,p) = bF^{x,p), (34) 

for a, h in the spectrum of A, is a {distributional, if the spectrum of A is continuous) pure state of the 
form (31). The associated wave functions V'l (x) , ■02 (2^) satisfy the eigenvalue equations A'il)i{x) — 
aipi{x) and A''ijj2{x) = b*il!2{x). 

Proof: A generic linear operator can be displayed in a fully symmetrized form through the Weyl 
prescription: 

Mq,p) =Jd^ldv a{^,rj)e'^^+'^^, (35) 
where a{^, rj) = ^Tr [A {q,p) e-'^^-'ip}. The Weyl transform of A is given by (cf.(12)): 

A{x,p) = Jd^Jdri a{^,ri)e'^''+'''P. (36) 

We now assume that F^{x,p) satisfies the stargenvalue equation (34) with A{x,p) given by (36). 
Using eqs. (14,29) we get for the first stargenvale equation (34): 

aF^{x,p) = ^ j di j d7] j dy a{i,7])e'^^+''^^^-'^^^h-'y^+'^^^^Z{x,y) 

= ^jd^jdrjjdy Q;(e,r?)e^«(^+^^/2) ^-vdy^-iyp e'^^^-Z{x,y) 

= ^jd^jdvjdy a{^,ri)e'^''-'yPe''^ye'^^y/^e-''^ye''^ye'>'^^'-Z{x,y), (37) 

where in the last step we integrated by parts in y and used the relation e~'^^ye^^y — 1. The 
Baker-Campbell-Hausdorff relation yields: 
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and substituting in (37) we get: 

We notice that the commutator of C = i^{x + hy/2) and D — r]h{^dx + ^dy) is a c-number and 
thus e'-'e^ = gC+Dgi[c,D]_ Substituting in (39) we get: 



e 'yPe 



Z{x,y). (40) 



Using (36) and (29) we finally obtain: 

h ,1 



A 



1 



Z{x,y) = aZ{x,y). 



Following exactly the same steps we get for the second equation in (34): 

1- 



Let us now define the new variables: 



u 



X + ^y 



Z{x,y) = bZ{x,y). 



d = + -d 



-d - -d 



yv = X- ^y 

and the new function G{u,v) = Z{^, ^) = Z{x,y). Prom eqs.(41,42) we have: 

A u, —ihdr 
A v,ihdy 



G{u,v) = aG{u,v) 
G{u,v) = bG{u,v). 



From where it follows that G{u,v) — (l){v)ip{u) where ip{u) satisfies: 



A 



u, —ihdu 



iIj{u) = aip{u), 



and (j){v) obeys to: 



V, ihdy 



(Piv) = b(f){v) ^ 



V, —ihdy 



(f)*{v) = b*(j)*{v). 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



Let ■ip2{x) = (t)*{x) and il)i{x) = ip{x). Then Z{x,y) = G{u,v) = ■?/'2(w)-?/'i(m) = 'ip2ix - ^y)ipi{x + ^y) 
and thus F^{x,p) is a (possibly distributional) pure state. Furthermore, t/'i and ip2 satisfy the 



proper eigenvalue equations. Finally if A is hermitian then eq.(46) reduces to A 
bip2{x).n 



X, —ihda 



1p2{x) 
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Later in this section we will provide the specification of the wave functions ■0i)'02 in terms of 
the original Wigncr function. First let us explore the physical meaning of the "initial conditions" 
for eq.(30). We define the mean n-th order momentum as: 



<p"(x,i) >12E 



Vl2{x,t) 



where 



Vi2{x,t) ^ jdp F^{x,p,t)^r2{x:t)tlj,{x,t). 



(47) 



(48) 



The function Z{z,j,t) then works as a generating function of the higher momentum correlations 

A'gp{x,ty. 



Vi2{x,t) = Z{x,0,t), 

i" dj 



A^I$p{x,t)^j:§,\nZ{x,j,t) 



(49) 



j=0 



In particular, we have: 



I A'^^^p{x,t) =<p{x,t) >12, 
\ Afjp{x,t) =<p'^{x,t) >i2 -<p{x,t) > 

From eq. (30,49) we can obtain the following recursive relation: 



2 

12 ■ 



(50) 



A 



(n+2) 
12 



p{x,t) 



'1 

2i 



A^;ip{x,t), Vn>l. 



(51) 



This means that we can compute all Ai2 p from A12P and Ai2P- Consequently, we get for n > 1: 

' AS?p(x,i) = (|)"£lnn2(x,i), if n is even 



n-l 



(52) 



where 



, , hi d ^ 

<p{x,t) >i2=y^ln 



^" r < p{x, t) >i2, if n is odd 

>^(a;,t)' 



(53) 



/tpi{x,t)_ 

A pure state could be seen as one where all the higher momentum correlations depend only upon 
A^-^^p and A^'^^p. This is not surprising since Vuix.t) and < p{x,t) >i2 provide suitable "initial 
conditions" for cq.(30) and thus completely determine the resulting pure state. 

We proved the consistency of the pure state quantum condition with both the Moyal and star- 
genvalue equations and studied the physical meaning of the "initial conditions of eq.(30). We now 
address the issue of deriving explicit expressions for the wave functions in terms of the original pure 
state Wigner function. This will be the purpose of the two following subsections. 
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4.1 Time independent pure states 

Let us assume that the time independent function F^{x,p) is square integrable and satisfies eq.(30). 
Then we already know that takes the form (31). How can we construct the functions '02 
from the knowledge of F^7 The answer is trivial: 

' ^^1(0:) = N,Jdp e^f^/^Fi^ (f ,p) = N^Z (f , f ) , 

(54) 

M^) =N2ldp e^f-/^ [f^ {I,p)Y = N2Z* (f, f) , 

where A^i, N2 are some complex constants (obviously A^i = ['?/'2(0)]~'^ -^2 = ['V^i(0)]~^)- Fur- 
thermore and since Z (^|, |^ is a square integrable function, then so are ipi, 1^2- There may be a 
problem if ipi{0) = or ip2{0) = 0. However, this can be easily overcome, as we now argue. Suppose 
that when we compute ^'1(2^) according to equation (54), we get zero identically. This is a sign that 
-02(0) = 0. In that case, we can use the alternative formula: 

M^) ^M.jdp e'^^'^^FY2 , (55) 

where Mi is another complex constant. Indeed, from (31), the right-hand side of (55) reads: 

^[V^^*(0)^i(a:) + ^;(0)^i(x)]. (56) 

In all cases of interest, the function ■02(3:^) and its derivative '02 (^) do not vanish simultaneously at 
the same point^. We then have ■02 (0) = 0, but ■02(0) 7^ 0; and hence (56) yields: 

(57) 

And thus the previous arguments still apply. These technicalities could be formally averted, if 
instead of evaluating ■0i or ^02 at x = 0, we had chosen some other point xq such that ^'i(x'o) 7^ 
and ■02(2^0) 7^ 0. However it turns out that in most cases, the calculation is easier to perform at 
X = (see section 7). 

In theorem 2 we proved that if F^ satisfies the generic *-genvalue equations (34) then F^ is a 
pure state and the associated wave functions satisfy the corresponding eigenvalue equations. Now 
that we have explicit formulae for the associated wave functions let us show that they provide the 
proper eigenfunctions for the case of the energy *-genvalue equation: 

Theorem 3: If F^ obeys the *-genvalue equations (27), then the associated wave functions -01, -02 
(eq.(54)) satisfy the corresponding energy eigenvalue equations (eq.(26)). 

Proof: From the first equation in (27) and from (14), we get: 

^For a solution of the Schrodinger equation, this would imply that '02 {x) vanish identically. 
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where V'^'^^x) is the n-th derivative of V{x). If we perform the substitution x — > |, multiply the 
equation by Nie^^^^^ and integrate over p, we get: Ai + A2 — A3, where: 



I dp £ [e^P^/^F^ (f,p)] = 



(59) 



2m 



= TV, / dp eW^' ^ (^)" (f )" (f ) Fi'T (f ,P) = V{x)M^), 



(60) 



^(x) = 7Vi£;i I dp e'P^/'^F^ (^,p) = E^Mx). (61) 
Assembling all the results, we get: 

-^^M^) + V{x)ij^{x) = ^iV^i(x). (62) 
Similarly, from the second equation in (27), we would get the following eigenvalue equation: 

^ '^\2{x) + V{x)iP2{x)^E2M^)- (63) 



2m dx^ 



It is easy to prove that the same results would hold even if we had used the alternative formula 
(55) in the case ■?/'2(0) = O.n 

This completes the time independent case. 

4.2 Time dependent pure states 

The time dependent case is shghtly more dehcate. If we apply naively the formula (54) to a time 
dependent Wigner function (31), we will get: ^^2(0) O'0i(^) t)- So there is an overall time dependent 
constant which cannot be distinguished from the own time dependence of the wavcfunction i/Ji{x, t). 
This is not surprising as the following argument illustrates. Consider a pair of states \'ipi(t) >, 
1^2 (^) >• From these we can construct the entries of the non diagonal density matrix element 
\ipi{t) >< ip2{t)\ in the position representation: 

p{x,y)^<x\ipi{t)><ip2{t)\y> . (64) 

Now suppose we multiply |'^i(t) > and \'4'2{t) > by the same time dependent phase eK"(*), then 
the density matrix element (64) remains unchanged. This means that, given some density matrix 
p, there are several wavefunctions (with different time evolutions) which yield the same physical 
predictions as does p. From among these we want to obtain the wavefunctions ipi, ip2 (unique up 
to complex multiplicative constants) which obey the Schrodinger equation. 
The desired functions are given by the following theorem. 
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Theorem 4: Suppose that F^{x,p,t) is a solution of the Moyal equation and of the pure state 
constraint (30) and let 



' i)i{x,t)^e-x^{\hiVi2{x,t) + {Jq dx' <p{x',t) >i2 

_ V'2(x, t) = exp {I \nVUx, t) + i Jo" dx' < p{x\ t) -ia*{t)} , 



where 



<p(0,t') > 



2m 



ait)^ I dt' 
Jo 

and Qi2 is the non-diagonal quantum potential: 

1 



^ + V{0) + Qu{0,t') 



Qi2(o,t) 



4mPi2 



dx^ 2V 



1 f&P 



12 



12 



9x 



(65) 



(66) 



(67) 



Then F^{x,p,t) can be expressed in terms of ipi and ■02 as the bilinear functional (31) and the 
functions ipi, ■02 obey the Schrodinger equation. 

Proof: If is a solution of the Moyal equation, 



dFl 

dt 



H,F, 



M 



p dF^ (-1)" d^''+^vd^''+^F^ 



(68) 



it is then easy to conclude that the dynamics of V12 and < p >i2 are dictated by: 

0, 



dPi2 I ^ 



Since Fj*^ is a pure state which satisfies eq.(30), then we also have (cf.(52)): 



(69) 



Afjp{x,t) ^ -'-^-^\nVu{x,t) 



It then follows that: 



d<p>i2 — 9 I h'^ 
dt dx 1 4m 



1 8^Pi2 _ 1 ( 1 dVi: 
ri2 dx^ 2 \Vi2 dx 



<P>1 
2m 



x 



= -|:{Qi2(x,t) 



<p{x,t)>\ 



2m 



+ V{x)] 



(70) 



(71) 



From these expressions we get: 

^,{x,t) = + i/^dx' <p(x',i) >i2 -ia(i)}V'i(x,i) = 



\ 2mJ 



dVi2 

2mVi2 dx 



^ r, 1 d<p>i2 _ i_ 

< P -^12 2m dx 



Qi2(x,t) + ^^^^ + n^) 



+ 



+i [Qi2(o, t) + + y (0)] - ^ [Qi2(0, t) + ^^M^ + V(0)] } t) 



(72) 



<p>12 dVi2 1 '9<P>12 _ i 

2mPi2 dx 2m dx h 



'Ql2{x,t) 



<p(x,t)>l 



2m 



■ + V(x)]}Vi(x,t). 
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On the other hand we have: 



1 &P12 



2Vi2 dx 



+ I < P >12 



+ 



1 a^-Pi2 

2Vi2 dx^ 



(73) 



<P>?2-tQl2(^,t) + i 



h V12 dx 



Comparing eqs.(72) and (73), we recover the Schrodinger equation (23). A similar result can be 

obtained for ip2{x, t). 

It remains to show that can be expressed in terms of the functions ipi, ip2 (65) according to 
(31). Since is a pure state we already know that: 



F^' 



12 



x,p,t) 



1 

2^ 



where A and B are some complex functions. Substituting in eqs.(65), we get: 

' M^,t) = i?(a:,t)exp{iln - fa(t)} 

< 

_ r2{x,t) = A{x,t)exp{-'^ln[^]+^a{t)} 



We then conclude that: 



hy 



hy 



hy 



hy 



ipUx -,t\ipi\x+ —,t \ = Alx -,t\ B [x + —,t 



(74) 



(75) 



(76) 



and it follows that F^ is indeed the bilinear functional of ipi and ip2 (31). □ 

Notice that the equations (65) are equally valid for the time independent case. However, the 

formula (54) is usually easier to evaluate, and it does pay off to use it in the time independent case. 
Before we end this section let us make some remarks: 



1) Equations (30) and (54) were introduced by Tatarskii in ref.[12]. He admitted only diagonal 
state solutions to equation (30). Here we accept non-diagonal pure states as well, because these 
will be crucial for the construction of mixed states in section 5. Tatarskii considered (54) as the 
right formula for time dependent states. Here we consider the alternative formula (65) for the time 
dependent case, because otherwise the wavcfunction will have an undefined time dependent phase. 
Furthermore, we considered a general stargenvalue equation (34) and proved that its solutions are 
pure states associated to wave functions satisfying the proper eigenvalue equations. The explicit 
form of these wave functions is given by (54). We checked this explicitly for the energy stargen- 
functions (Theorem 3). In summary, the former results provide the complete generalization and 
specification of Baker's converse construction for pure state Wigner functions. 

Mathematically this specification can be implemented as follows (to make it simpler we consider 
only real pure states, the extension to the more general case being straightforward): let us denote 
by J-'^re the space of all real, square intcgrable and normalized functions defined on phase space 
and which satisfy the condition (30). This space is equivalent to the quantum Hilbert space of 
distinct physical states. Let [ip] be the equivalence class of normalized wave functions belonging 
to the Hilbert space Ti and such that (p G if and only if there is a global phase factor e'^ such 
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that = e*^V^. Let [H] be the set of these equivalence classes. Then there is a one to one map 
: [n] — > T^,, defined as: 

and displaying the following inverse map: (ji^^ : J^^^e — ^ [^1' with: 

F^{x,p) ^ (A^)~' {f^{x,p)) = [ip{x)] : V'(a;) ^ N j dp e^^^/^F^ , (78) 

where iV is a normalization constant (defined up to a global phase factor). Furthermore, if ip{x) 
is an eigenfunction then F^{x,p) is the corresponding stargenfunction and vice- versa. Finally, 
if ip{x,t) satisfies the Schrodinger equation then {[ip{x,t)]) satisfies the Moyal equation. Con- 
versely, if F^ {x,p,t) satisfies the Moyal equation then there is a wave function (f){x,t) G [■?/'(a;, t)] in 
the equivalence class of ip{x,t) (eq.(78)) that satisfies the Schrodinger equation. The explicit form 
of (j){x,t) is given by eq.(65). 

2) The pure state quantum condition (30) can be rederived in a different (but equivalent) functional 

form by using the momentum representation of the original wave functions as a starting point. For 
a generic quantum state {ip > the momentum i/j{p) =< p\i/j > and the position ijj{x) =< x\ip > 
representations are related by: 

i>{p) = f dx ij{x)e-'''P/^ and 7jj{x) = f dp Vi(p)e^^^'/^ (79) 

from where it follows an alternative functional form of the pure state Wigner function: 

FZrei^,p) = ^Jdy e-^y^r - y) ^ + y) = ^ / e-'^^'/^* {p-k)^{p + k). 

(80) 

We now follow the procedure of eqs. (29,30) and define: 

E{y,p) = jdx e'^yF'^{x,p), (81) 

and from (80) we easily conclude that a real, square integrable and normalized phase space function 
is a pure state Wigner function iff it satisfies: 

lnE(7/,p)= - — lnE(y,p). (82) 



dy'^ ' \2 ) dp 

Hence, the former equation provides an alternative (equivalent) formulation of the pure state quan- 
tum condition (30) that in some specific cases might be easier to evaluate than formula (30). Finally, 
we should also notice that the pure state conditions (29,30) and (81,82) can be rederived for arbi- 
trary non-commuting phase space observables A, B satisfying [A,B]m — 1- 

3) The pure state quantum condition implies (but it is not equivalent to) Heisenberg's uncertainty 
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principle. In the next section we will prove this result for the more general set of mixed quantum 
states (which includes as a subset the pure quantum states). Here, we follow Tatarskii through the 
following elucidative example of the relation between Heisenberg's uncertainty principle and the 
pure state quantum condition. Consider a Hamiltonian quadratic in the positions and momenta 
hke for instance the simple harmonic oscillator. The Moyal equation reduces to the classical Liou- 
ville equation and all references to Planck's constant thus disappear. Indeed the time independent 
Moyal equation reads: 

(83) 

Any function of H solves eq.(83). In particular, let us choose the following solution: 

F'^ix^p) = ^e-"^^^'^'), (84) 

ZTT 

where a is a positive real constant. The function F'^{x^p) is real, normalized and square integrable. 
However if we compute the position and momentum dispersions, we get: 



1 



{^2xY = 3::-:^, {^2pY = -, (85) 

which means that: 



Aax • A2P = — . (86) 
auj 

Since the constant a is arbitrary, there is no lower bound for the product in eq.(86), in which case 
Heisenberg's uncertainty principle may be violated. So (84) will in general not be an admissible 
quantum state. 

If, however, we impose the pure state condition (30), we get: 

and the product (86) will now respect Heisenberg's bound. This example also shows that the pure 

state quantum condition is more restrictive than the Hciscnberg uncertainty relations. In fact the 
Heisenberg bound implies only that a < We will came back to this issue in the next section. 

4) An interesting question is whether classical dynamics preserves the pure state quantum con- 
dition. We will now prove that in general this is not the case. To do so we consider the example of 
an anharmonic oscillator described by the Hamiltonian: 

ii"(x,p) = — + -mwV + Ax^ (88) 
2m 2 

where A is a positive and real coupling constant. The Moyal equation reads: 

^ = [H, F^U - -^T^ + {mw^- + 4Aa;^) ^ - Af^^^, (89) 
ot m ox ^ ^ op op-^ 



17 



and up to first order in t we have: 

F^{x,p, t) = F^ix,p) + t (^-^^ + (rnw-:, + 4Ax^) ^ - Xh'x^j + 0{t% (90) 

where the subscript Q stands for the quantum time evolution of F^{x^p). If we use the Poisson 
bracket instead and calculate the classical (Liouville) time evolution of we easily get: 

F^{x,p,t) = F^{x,p,t)+tXh'x-^ + 0{f). (91) 

To proceed we calculate Z{x,j,t) up to first order in t both for the classical and quantum time 
evolutions. We have: 



ZQ{x,j,t) = Jdp e^^PF^{x,p,t) = Zo{x,j)+itZ?{x,j) + 0{t') 
Zc{x,j,t)^Jdp e'^PF^{x,p,t) = Zo(x,j)+itZ^{x,j) + 0{t''), 



(92) 



where: 



Z?{x,j) = - (mw'xj + 4AxM Zo - h^XxfZo = Z^{x,j) - h^Xxf Zo{x, j), (93) 



and Zo = Z^ix^j) ^ J dp e'^PF^{x,p). Now notice that: 

t + 0{t^)^ In Zo{x J) +0{t'), (94) 



d 

In Zq{x, j, t) = In Zo{x, j) + ^ In Zq{x, j, t) 



<=o ^o(a;,i) 
and thus, up to the first order in t, the pure state quantum condition is written: 

That this equation is satisfied is a consequence of Theorem 1. This result can also be checked 
exphcitly using (93). Here we want to see if eq.(95) is also satisfied if we consider the classical 
evolution instead. We have: 

^ ( z^{x,j) \ ^(ny^ ( zr{x,j) \ 

df[Zo{x,j)) [2) dx^[Zo{x,j))^ ^ ^ 

6h^Xxj = 0, 



92 (z?{x,j) + h'XxfZoix,j)\ (hV {z?ix,j) + h'XxfZo{x,j)\ _ 



df \ Zo{x,j) J \2j dx^ V Zo{x,j) J 

where we used eq.(93). The last identity is satisfied only if A = in which case the anharmonic os- 
cillator reduces to the simple harmonic oscillator for which the Moyal and Liouville time evolutions 
coincide. We conclude that in general the classical dynamical structure is not compatible with the 
pure state quantum condition. An interesting question for future research is whether the quantum 
condition uniquely determines the Moyal dynamics. 

5) Another interesting result is that the classical measurement procedure does not preserve the 
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pure state quantum condition either. Assume that F^{x,p) is a classical phase space distribution 
(i.e. is real, positive defined and normalized). There are in fact few cases where these conditions 
are compatible with the phase space quantum condition (the problem, as is well known, is posi- 
tivity). The Gaussian distribution eq. (84,87) is one of those few examples that can be regarded as 
both a quantum and classical distribution. So assume for the moment that F^{x,p) is a Gaussian 
distribution. Under a classical measurement of a generic classical observable A{x,p) with output a, 
F^{x,p) will collapse to: 



where is a normalization constant. The new distribution is still well defined as a classical distribu- 
tion (is positive, real and normalized) and displays a well known feature of classical measurements: 
it preserves the original functional form of F^{x,p) on the phase space hypersurface A{x,p) — a. 
This feature is of course not shared by the quantum measurement and it is trivial to realize that 
in general the outcoming distribution (97) is not a pure (even distributional) quantum state of the 
form (17) (take for instance A{x,p) = p). In fact, the distribution (97) in general does not satisfy 
the Heisenberg uncertainty relations (for A{x,p) = p the Heisenberg uncertainty relations would 
imply that F^{x, a) — Const, \/x) and thus cannot be a pure (distributional) quantum state. We 
conclude that the preservation of the pure state quantum conditions through the measurement pro- 
cedure imposes a reformulation of the classical collapse prescription. Again, an interesting question 
is whether the quantum conditions uniquely determine the standard form of the quantum collapse 



In this section we generalize the previous results to mixed states. In particular we want to deter- 
mine when does a phase space function F^{x,p) represent an admissible mixed quantum state. As 
always it must be real and normalized. But now we want to find some suitable condition which 
plays the role of (30). 

Definition: Let F^^^{x,p) be an arbitrary pure state Wigner function, i.e. a real, square in- 
tegrable and normalized function which obeys eq.(30). Then F^{x,p) is an admissible quantum 
Wigner function if it is real, normalized and satisfies the following " mixed state quantum condition" : 



for ah pure state Wigner functions F^j.^{x,p) G ^^r-e-a 
We shall now prove the following Lemma: 

Lemma 3: A real and normalized phase space function F^[x,p) is of the form (21) iff it obeys the 

"mixed state quantum condition" (98). 

Proof: Suppose that the real and normalized function F^{x,p) satisfies the inequality (98). Let 
us define the candidate density matrix in the position representation: 




(97) 



[31]. 



5 Mixed states 




(98) 




(99) 
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where (at this point) y) might only make sense as a distribution. We have: 

Trp ^ J dx p{x,x) ^ J dx J dp F^{x,p) = 1, (100) 

p{x,yr ^ jdp e-TMF^ (^'^^) = Ml/,^), (101) 

and so p is hermitian and normahzed. Moreover, for an arbitrary normahzed wavefunction 'il){x), 
we have: 

Jdujdv ij*{u)p{u,v)^{v) = Jdpjdufdv e-f ("-"^F^ ^*(M)^(t;) = 

^hjdpjdxjdy e-'yPF^{x,p)r {x - ^) {x + 'f) ^ 2nh J dx J dp {x,p)F^,,{x,p) > 0, 

(102) 

where we used (98) and: 

FZreix,p) -^Jdy e-^y^r - y) + y) ■ (103) 

We also performed the substitution: 

u — x — , V — X+—. (lU4j 

Since ijj{x) is arbitrary, wc conclude that the quadratic form associated with the matrix p is non 
negative. It follows that p satisfies all the postulates of the quantum density matrix discussed in 
section 2. Consequently, it can be expressed as a superposition of pure states as in (10) and it 
satisfies Trpp < 1 =^ J dxdy p*{x,y)p{x,y) < 1. Hence, we can invert the Fourier transform (99) 
and conclude that F^{x,p) is of the form (21). 

Conversely, if F^{x,p) satisfies eq.(21) then F^{x,p) is real, normalized and satisfies: 

dx dp F^{x,p)F^^^{x,p) = ^pf) dx dp F^^^^p{x,p)F^^^{x,p). (105) 
J J J J 

Let ifj ^{x) , (t){x) be the wave functions associated to F^^^ f^{x,p) and F^^^{x,p), respectively Both 
Fpure,i3i^^P) Fpurei^^P) display the standard functional form given by equation (103) and we 
have: 

Jdxjdp F^^^^^ix, p)F^,,{x, p) = ^jdxjdyrp{x- f ) (x + f ) 0* - f ) (x + f ) = 
= ^Jdujdv ilj}){u)iljf3{v)(p* {v)(p{u) = 2^ \jdu ip*^ (u) (p {u)\^ > 0, 

(106) 

where we used eq.(104). Substituting this result in eq.(105) we conclude that F^{x,p) satisfies the 
mixed state quantum condition (98). □ 

Corollary: All pure quantum states obey the mixed state quantum condition. 
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Proof: This result follows immediately from eq.(106).n 



Theorem 5: The mixed state quantum condition is preserved through the Moyal time evolution. 

Proof: We want to prove that if condition (98) is verified ai t — 0, then it will hold at all later 
times, provided F^(x,p) evolves in accordance with Moyal's equation. This a trivial result that 
follows from theorem 1 and the linearity of the Moyal equation applied to a generic mixed state of 
the form (21). Alternatively, we shall provide a more direct proof that explicitly exhibits the time 
evolution of the mixed state quantum condition. Let then: 

C/(i) = e^^^^^^'^^ (107) 

be the Weyl map of the quantum evolution operator. The noncommutative exponential is defined 
by: 

(108) 

no{x,p) = i, nn+i{x,p) = A{x,p) *nn{x,p), n>o. 

Then the solution of Moyal's equation (24) is: 

F^{x,p, t) = U(t) * F^(x,p, 0) * U{-t). (109) 

Using the identity [11] 

J dx J dp A{x,p) *B{x,p) = J dx j dp A{x,p)B{x,p), (HO) 

we conclude that: 

Jdxjdp F'^{x,p,t)F^,,{x,p)^JdxJdp U{t) * F^{x,p,0) *U{-t) * F^^^{x,p) ^ 
= JdxJdp F^{x,p, 0) * U{-t) * F^,^{x,p) * U{t) = Jdxjdp F'^{x,p, Q)F^rei^,p, -t). 

(Ill) 

Since the evolution of a pure state is a pure state, from (98) follows that: 

Jdxjdp F'^{x,p,t)F^,,{x,p) > 0, Vt e 5?, (112) 

and thus we conclude that F^{x,p, t) will be of the form (21) at all times. □ 

The previous analysis leads to the explicit form of the one to one maps relating the space of 
physical density matrices acting on the quantum Hilbert space and the space of mixed state Wigner 
functions. These maps provide a suitable generalization of the ones previously presented in the 
context of pure state Wigner functions eqs. (77,78). Let us denote by J^^ixed the space of real and 
normalized functions in phase space which obey the mixed state quantum condition (98). Likewise 
let Dili.) be the space of linear operators acting on the Hilbert space Ti which obey to the conditions 
1) to 3) of section 2. The two spaces are related by the Weyl map: Vw : 'D{H) — > ^mixed^ which 
is given by: 

m ^ im) = F^i^,p, t)-^Jdy < ^ + Y^mi^ " ^ ^^^^^ 
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and is one to one. The inverse map is : J^mixed — ^ ^(^)) reads: 

F^{x,p,t) ^ V^^ (F^{x,p,t)) = p{t) = n j dx j dy j dp e'yPF^{x,p,t)\x 



Obviously, we have J^^^^ C J^Zxed- 



(114) 



pure 

Now that we have completely specified the relation between the kinematical structures of the 
Wigner and the density matrix formulations and proved their formal equivalence, we can import 
the notion of complitude from the space V{T-C) to the quantum phase space J-'Zixed subsequently 
address the related issue of which quasi-distributions appear with which weight in a given mixed 
state Wigner function (21). 

Suppose that {|a >, a e /} is a set of orthonormal eigenstates of some operator A which is 
complete (i.e it generates the full Hilbcrt space Ti). All admissible quantum density matrices 
are thus expandable in the basis {\a >< b\, a,b & I}. These basis elements satisfy the following 
orthogonality conditions: 

\a >< b\c >< d\ ^ 5b^c\a >< d\, \/a,b,c,deI. (115) 
If we multiply this equation by j^^^ and apply Weyl's map, we obtain: 

FTb*FZ = ^5,,cFZ, ya,b,c,del, (116) 

where F^(x,p) are the non diagonal Wigner functions -^Vw {\a >< b\). They obey the *-genvalue 
equations: 

Aix,p) * FZix,p) = aFZix,p), FZ{x,p) * Aix,p) = &F„^(x,p), (117) 

where A{x,p) = Vw Because of the one to one correspondence between V{H) and J-'mixed^ we 

conclude that |F^(a;,p), a,6 G /j generates J^Zxed- This means that for any F^ e ^mixed-> "^^ 
have: 

F^(x,p) = ^ PabFZ{x,p), Pab e C. (118) 

a,b&I 

From (116,118), we have: 

FZ^ * F^ * FZ = Y:jab5m,a5b,rFZ = (^P-^nT- (US) 

And so: 

p^, = {2nnf - ^. (120) 

ns 

A trivial calculation shows that these are just the elements of the density matrix in the basis 
{\a >, a e I}: 



< 



r\p\m >— J dx J dy < r\x >< y\m > p{x,y) — J dx J dy J dp %l)*{x)'il)„i{y)^ '^^^ y)p^ ^2±a^p^ - 

= ^ Y.a,ba Pab Idxjdyjdzjdp rr{x)^m{y)e-'^^'-'^+''^ra {^^) ^b = 

X^a,b6/ Pab^r,b^m,a Pmr i 

(121) 
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where ipnix) —< x\n >. 

Alternatively we may expand F^{x,p) in the form F^{x,p) = J2fS(^BPi3F^re,i3i^TP) where the 
associated pure states are the solutions of the left and right stargenvalue equations: 



This result follows immediately from applying the Weyl map to the diagonal representation of p 
(eq.(lO)). 

To proceed we provide an alternative form of the mixed state quantum condition. 

Theorem 6: The mixed state quantum condition is equivalent to the condition: 

< 9*{x:P)*g{x,p) >FW= j dx j dp {g*{x,p) * g{x,p))F^ {x,p) > 0, (122) 

to be satisfied for all complex valued functions g{x,p), provided the integral exists. 

Proof: We begin by proving that if F^{x,p) is a pure state then it satisfies eq.(122). We have: 

j dx J dp {g*{x,p)*g{x,p))F^^^{x,p) 
= 27rh jdxjdp {g*{x,p) * g{x,p)) * * F^,^{x,p)) 

= 27rh J dx j dp {g{x,p) * F^^^{x,p)y * {g{x,p) * F^^^{x,p)) 
= Jdxjdp \g{x,p) * > 0, (123) 

where we use eq.(llO) and the identity (28) for pure states. The generalization of the former result 
for mixed states is straightforward. Using eq.(21) we have: 

[dx [dp ig*ix,p)*g{x,p))F'^ix,p) = pJ dx [ dp ig*ix,p)*gix,p))F^,,^p{x,p) > 0, (124) 

where we used the result (123) for each of the pure states F^^^^p{x,p). 
To prove the converse result we use eq.(28) and get: 

Jdxjdp F^^^{x,p)F'^{x,p) = 2nh Jdxjdp {[F^rei^,p)]* * F^^^{x,p)) F'^{x,p) > 0, (125) 

where we also used the reality of the pure state Wigner function and the condition (122) for g{x,p) = 

The condition (122) was introduced in a different context by C. Zachos and T. Curtright and 
proved to imply the Heisenberg uncertainty relations for arbitrary non-commuting observables [35]. 
Hence, we conclude that the mixed state quantum condition implies the Heisenberg uncertainty 
relations. Finally we prove that: 

Theorem 7: Let Q be the set of complex value phase space functions g{x,p) that satisfy the 
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normalization condition: / dxdp g*{x,p)g{x,p) — 1. Then the set of mixed state Wigner functions 
is given by: 

^Z.ed = {g*{x,p)*g{x,p) : g{x,p) e S}. (126) 



Proof: We begin by proving that all phase space functions of the form g*{x,p) * g{x,p) with 
g{x,p) G Q are mixed quantum states. In fact g*{x,p) *g{x,p) is real, normalized and by virtue of 
theorem 6 it satisfies: 

Jdxjdp F'^ix,p) {g*ix,p) * g{x,p)} > 0, VF'^ G J'Z.ed ^ 

Jdxjdp F^,^{x,p) {g*{x,p) * g{x,p)} > 0, ^F^^^ G ^^t,. (127) 

From Lemma 3 we conclude that g*{x,p) * g{x,p) is a mixed state. 

To prove the converse result let us start by recalling that an arbitrary mixed state can be cast 
in the form: 

where the pure states satisfy: F^^^ f^{x,p) * F^^^ ,^{x,p) = Fp^^g_^(x,p). Let then: 



9{^,P) = E ^/^^^PpF^reA^^P)^ (129) 
peB 

and we have: 

g*{x,p)*g{x,p) = \/'^^^P»\/'^^^PpFZre,ai^^P) * Fpured^^P) (130) 

= E VP^^^,pF^reA^^P) = E PpF^reA^^P) = FZi^ed{x,p). 
a,/3eB f3eB 

We conclude that an arbitrary mixed state Wigner function (128) can be written as in (126). Finally, 
notice that the decomposition g*{x,p) * g{x,p) is defined up to an unitary transformation, i.e. if 
g{x,p) satisfies (130) then U{x,p) * g{x,p) also satisfies (130) for generic unitary symbol U{x,p).a 

Equation (126) provides an elegant, "covariant" characterization of the space of physical states 
in quantum phase space. An interesting feature of this equation is that it only depends on the 
starproduct defined on the phase space. This leads to the conjecture that this characterization 
(126) might also be valid in the more general setting of curved phase space manifolds [31]. 

Some remarks are now in order: 

(i) First, notice that the mixed state quantum condition, besides implying the Heisenberg uncer- 
tainty relations, also implies the positivity of the marginal distributions. In particular we have: 

/ dp F'^{x,p) ^J2pp [dp F^reA^,p) = E > 0, Vx G 3?, (131) 

PeB ■' peB 

which (together with the normalization condition) means that we can interpret the marginal dis- 
tribution as a true probability density for x. Likewise, we can show that: 



/ 



dx F^{x,p) > 0, Vp G 3?, (132) 
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and thus interpret this quantity as the probabihty density for p. 

(ii) The properties of Wigner functions associated with quasi-distributions have been analyzed in 
the past. In particular, it was shown [12] that a real and normalized pure state satisfies: 



j dx j dp 



^P^e(^,P)]' = 2^- (133) 



Likewise, for mixed states we have: 



J dx j dp 



FZU^,p)]' < (134) 



Although these conditions can certainly be regarded as preliminary tests to eliminate spurious 
states, they are by no means sufficient conditions. The necessary and sufficient conditions for real 

normalized phase space functions to be pure or mixed quantum states are those stated in equations 
(30,98). Conditions (133,134) are a corollary of the latter. For instance eq.(133) is a straightforward 
consequence of eq.(28). 

(iii) Another interesting fact is the following. We want to show that if is a pure state and 
solution of the time-independent Moyal equation, then it must be of the form: 

F^{x,p)^FZ{x,p), (135) 

for some n, where |f^(x,p), m G /| are the diagonal *-genstates of the energy: 

^ * ^1. = ^In * ^ = EmFZm- (136) 

Let then F^ — J2n mei PnmFZn- If it ^ solution of the time-independent Moyal equation, then 
we have: 

Y: Pnm{En-Er,)FZ,^0. (137) 

Since ^F^{x,p), n, m e /j is a complete basis, we conclude that 

F^(x,p) = Yl PmmFZm{x,p). (138) 

If in addition we impose that F^ be a pure state, then = except for just one m E I, and we 
recover (135). This is exactly, what happened in our example in section 4. By imposing the pure 
state condition, we forced F^ to be a diagonal *-genstate (84,87) (in this case the fundamental 
state of the harmonic oscillator). 

(iv) Finally, it is instructive to return to the example discussed in section 4. There we saw that 
by imposing the pure state quantum condition we obtained a state that satisfied the Heisenberg 
relations. In fact these were saturated. Now we may admit states of the form (84) that are mixed. 
In that case we have to impose the mixed state quantum condition (98). 

Let us then suppose that the space ^Zxed is generated by the *-genstates of the simple harmonic 
oscillator. If F^ given by eq.(84) is a mixed state then, in particular, it must satisfy the inequality: 

J dx j dp F^{x,p)F^{x,p) > 0, (139) 
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where 

The diagonal Wigner function is the first excited state of the simple harmonic oscillator. It is 
a solution of the *-genvalue equations: 

H{x,p) * F^{x,p) = F^{x,p) * H{x,p) = E^F^,{x,p), (141) 

where Ei — ^uU. Upon substitution of (140) in (139) we get: 

/d./dp = (A - a) > 0. (142) 

We conclude that: 




(143) 



and the product (86) automatically satisfies the Heisenberg relations. 

In the case of a pure state the quantum nature of the state is made explicit through the appear- 
ance of Planck's constant in the constraint (30). It seems that the mixed state quantum condition 
(98) is not as explicit. However, if we recall that the pure states appearing in (98) must satisfy 
(30), then the quantum nature of the mixed state is no longer elusive. 



6 Alternative quasi-distributions 

The Wigner function is the quasi-distribution associated with the Wcyl correspondence rule. This 
means the following. We can use the Wigner function to compute e.g. the mean value of a given 
operator A (eq.(18)) if we choose to cast A in a fully symmetric form (Weyl order). Alternatively, 
we may decide to impose that the momenta precede the positions (standard order) or vice-versa 
(anti-standard order). But in that case, in order to keep the mean value unchanged, we must alter 
the quasi-distribution. Consider for instance the following exponential: e*^'^"'"*''^. If we expand it, 
all powers will be completely symmetrized in q and p. We conclude that the exponential is written 
in the Weyl order. Using the Baker-Campbell-Hausdorff formula we may express it in a different 
form: e^^^e^'^e'^^^ . Since the positions precede the momenta, the exponential is now written in 
the standard order. On the contrary, we may choose to place the momenta before the positions 
(anti-standard order), in which case we would get: e*''^e*^'^e~~^''. 

The previous analysis can be generalized in the following way. The exponential can be written 
in an arbitrary order according to: 

(144) 

where f{C,v) is called the Cohen function [1, 36] associated with the ordering prescription in ques- 
tion. We assume it to be an analytic function for all {C,r]). For instance, in the previous examples 
the Cohen function amounts to fs{C,v) — eT^^'^ (standard order), fAsiC^v) — ^^^^ (anti-standard 
order) and fw{^,v) = 1 (Weyl order). 

For a given correspondence rule (i.e. f{i,ri)), the quasi-distribution reads: 

Ff(x,p,t) -^^Jd^JdriJdx' < x' + f\m\^' - y > mv)e'^^''-'^-''''- (145) 
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We are now faced with the following problem. How can we tell whether a given function belonging 

to the space of all functions living in a phase-space with non-commutative product and Cohen's 
function f{C,,Ti), corresponds to a physical quantum state, i.e. belongs to J^^^ixed^ Moreover, if it is 
a physical state, is it pure or mixed? And, of course, we would like to reconstruct the wavefunction 
or the density matrix from the quasi-distribution. A generic quasi-distribution may have very 
different properties from those of the Wigner function. For instance, the Metha function F^{x,p, t), 
associated with the standard order, is not a real function. 

To get a feeling of the kind of difficulties one faces, let us consider the "mixed state quantum 
condition" (98). If we choose the standard order we would write this condition as: 

Jdxjdp F\x,p)F^l,{x,p) > 0, (146) 

where F^^^^^ is the quasi-distribution of an arbitrary pure state in the anti-standard order. So if we 
are interested in analyzing the properties of states in J-'mixed^ forced to study J-'^^^e well. 

The Wigner function corresponds to the simplest case because fwi^,v) = fwi^^v) = 1- 

The shortest route to solving this problem is the following. Prom eq.(145), we have the following 
relation between two quasi-distributions with Cohen's functions fi and /2: 

F\x,p,t) = ^ldx'l dp' I d^ I drj e^«(-'-)+^''(^''-f)||^F2(x',p',t). (147) 
In particular, if we choose f2 — f and /i = fw = 1, we get: 

F'^{x,p,t) = ^/ dx'l dp' I dC I drj e'^^^'-^^+'^^'-P^r\C,v)Ff{x',p',t). (148) 
The momentum generating function (29) then reads: 

Z{x,j,t) ^ jdp e^^PF'^{x,p,t) = ^J dx'l dp' J di e'^^^'-'^+'^P'r\C,j)F^{x',p',t). (149) 

Having calculated F^ and Z the whole analysis described in sections 4 and 5 can be directly 
applied. All the formulae presented depend only on Z and F^ . We can then determine whether 
Ff is a physical state, whether it is pure or mixed and we can construct the wavefunction or the 
density matrix from it. 

In the remainder of this section we will just state the explicit formulae for the most popular 
quasi-distributions found in the literature. 

For the Metha function {fs{C,,'r]) = e~^^^), we have: 



(150) 



F^{x,p,t) ^^Jdx' Jdp' e-W-^W-v)F^{x',p',t), 
Z{x,j,t) = Jdp e-'^PF^ (x+f,p,t). 
Similarly, for the Kirkwood function (^fAs(C,v) — e^^^), and Z read: 

' F^{x,p,t) = ^Jdx' Jdp' eir(^'-^)(f'-f)F^^(x',p',t), 

Z{x, J, t) = Jdp e'^PF^^ (x+f,p,t). 



(151) 
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The P-function of Glauber and Sudarshan is associated with normal ordering. Let us define the 
creation and annihilation operators of the simple harmonic oscillator: 

a — . {mujq + ip) , = . {mujq — ip) . (152) 

Normal ordering consists in writing the creation operators a} before the annihilation operators a. 
The associated Cohen function is fN{^,v) = e*^~^ 4^. We then get: 

(153) 

Z{x,j, t) = ^e-^j' / dx' J dp' e-'T^^'-^^'+^^f'F^ {x',p', t) . 

The Husimi function is defined by fH{^-v) = e~^^ 4^, where k is some positive real constant 
with dimensions of a frequency. It then follows: 

' F^{x,p,t)^^Jdx'Jdp'JdUdv e^+^^(-'-)e^+^''(^'-^')F^(x',p',i), 

(154) 

. Z{x,j,t)^^e'-^^'jdx'Jdp'!dC e^+^«(-'--Wiri^(x',p',i). 
In particular for n = u,we recover the Q-function F"^'^(x,p, t) associated with anti-normal ordering. 

7 Example 1: the simple harmonic oscillator 

To illustrate the results of section 4, let us consider the time independent harmonic oscillator. For 
simplicity, we choose units such that h — 1 and take m — cu — 1. The Hamiltonian reads: 

H(x,p)^^ + -. (155) 

In refs.[5, 9, 11] it was shown that the normalized diagonal *-genfunctions are given by: 

FZ{x,p) = 2(-l)"e-2^L„(4i/), n = 0, 1, 2, • • • , (156) 
where Ln{x) are the Laguerre polynomials. They are solutions of the *-genvalue equations: 

H{x,p) * FZ{x,p) = FZ{x,p) * H{x,p) = E^FZ{x,p), (157) 

with 

= n + ^, n = 0, 1, 2, • • • (158) 
It was also shown that the Laguerre polynomials admit the following integral representation [9]: 

L^^x)^^^^<f dz {z^ + l)-^z-^''-^exp(^^], n = 0,l,2,---, (159) 



27ri Jc \z^ + 1 
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where C is a closed contour such that 2; = is the only singularity in its interior. Substituting 
(156,159) in (54) we get: 



l^-2n-l 



X 



X J dp e*^^exp 
If n is odd we get {n — 2m + 1): 

V'2m+l(3^) = 



(160) 



2m+l 



^^/7^ 



e 2 (f) dz 



c Vl - z' 



exp 



(161) 



Applying Cauchy's residue theorem, we conclude that the contour integral in the previous equation 
vanishes. Consequently, the odd modes will have to be evaluated according to the alternative 
formula (55). 

Let us first compute the even terms (n — 2m): 



iP2m{x) = f dz exp 

l^J-K JC VI 



2 \l-z^ 



N2m -ai/-, / \ 



where 



4m— 1 



C2m{x) ^ f dz ^ exp 



x^z'^ 



(162) 



(163) 



If we expand C2m(3^) in powers of x using Cauchy's theorem, we get. 



27ri(-l)' 



22m^2m ^ Y^,^_^^k^2ra-k ^2{m-k) j 2m j . ^ . 3 . 5 . . . - 1) I . (164) 



fc=l 



2A; 



The term in curly brackets is readily identified as the Hermite polynomial of order 2m [37] : 



C2m(a;) = ^,22m-l -^2m(a:). 



(165) 



Consequently, equation (162) reads: 



V'2m(2;) = _yc^2-m-\ ^ ' -^2m(a;) = 



(-1)" 



m 



2"^ [0F(2m)!] 



1/2 ^2m('^) ) 



(166) 



where the constant N2m. was determined by requiring that -02^1 be normalized and hjyx) is the 
hermite function: hnix) = (— l)"e"'^^/^/J„(a;) = e^^/^^e""^^. 

It is interesting to note that, from this procedure, we obtain an integral representation for the 
even Hermite polynomials: 



H2m{x)^m\A^{-ir 



dz z-^'^-^ 



exp 



2 4 
X Z 



c 2m y/T^ ' \z^-lj' 



(167) 
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and also a relation between Laguerre and Hermite polynomials: 



Jdy e(^+'y^'L2m {2{x' + y'))^ 
Let us now try to compute the odd solutions: 



m!4" 



-H2m{'2x). 



(168) 



(169) 



-'-^fcdz (fei) (1 - zr'^' exp [-f (i±$) 

i^c^^ (1 - ^r'^" (1 + exp [-f (i^s); 

A simple calculation shows that the second integral in the previous expression vanishes. On the 
other hand we can easily show that: 



(i 2xN^. 
ax it: jc 



f^dz ^-4— 1(1-^4) '^'exp 



x^ l + z^ 



2 VI --^^ 



(170) 



Comparing (169) and (170), we conclude that: 

7r(-l)-+i7VL+i 



^^2m(2;) - Xh2m{x) 



(171) 



It is well known that the Hermite functions satisfy the equations [37]: /i^ = xhn + hn+i- Conse- 
quently, we get the following normalized solution: 



-1 



\m+l 



h2m+l{x)- 



(172) 



2"* [2V7r(2m + 1)!]^^^ 
Moreover, we also get an integral representation for the odd Hermite polynomials: 

H2rr.^^ix) = 2x4-m!(-l)- £ ^ z"^-^ (l - z')''^" exp . (173) 

From (167,173) we can construct an unified expression for the Hermite polynomials Hn{x) {n — 
0,1,2,- ■■): 



Hn{x) 



{2iT 



'2n+ (-l)'^- 1^ 



7 ( — 1)" 

Jc 2m ^ ' 



(174) 
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8 Example 2: the free Gaussian wave packet 



Let us consider the following real Gaussian Wigner function at time t = 



1 



F {x,p, 0) = ~^6xp 



X 



2al 



Now suppose that it evolves according to the free Moyal equation {V{x) — 0): 



dt 



2m' 



M 



p dF^ 
m dx 



The solution is: 



F^{x,p,t)^F'^ {x--t,p,0 



w 



P 



1 



m y nh 
A straightforward calculation shows that: 



exp 



V{x,t) = (27ra2)"^/^exp 



2at 



< p{x,t) >=Po + 



h?t{x—ut) 



Am 



1 {x~u€f 



2at 



where 



Po 
m 



1 + 



If we write in the polar form, ip = Re^^^^ (where S and R are real functions), 
(65,177,178): 



R{x,t) = [V{x,t)]^^^ = (27ra^) ^^%xp 



(x - uty 



and 



S{x, t) = dx' < p{x', t) > -a{t) 



n^tx 

















1 uH'^' 


^ 2m 


Po - 






1 2m 









Using the fact that = ^^2^^^^ , we get: 



olx, t) = — arctan 



1 \ n^t , 

2 I +Po\x - -ut ) + ^ g _g (a: - Mt)' 



0, 



Prom (180,182) we finally get: 



ip{x, t) — (^TTSf^ ^ exp 



(a; — uty i / 1 

+ j-Po [x - -ut 
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where 



St = (7o 1 + 



iht 



) 



(184) 



This is weU known to be the time evolution of the free Gaussian wave packet [38] . 

9 Conclusions 

Let us briefly summarize the results of this paper and comment on some possible future research: 
we determined necessary and sufficient conditions for pure and mixed phase space quantum states. 
We showed that this kinematical structure is fully compatible with Wigner quantum mechanics and 
used the emerging formalism to completely rederive standard operator quantum mechanics from 
the phase space formulation. These results extend and fully specify Baker's converse construction 
both to a general stargenvalue equation (associated to a generic phase space symbol) and to the 
time dependent case. 

Furthermore, we proved that the quantum conditions imply the Heisenberg uncertainty relations, 
a correction to the classical Liouville equation and to the classical measurement procedure. These 
results lead to the speculation of whether the quantum conditions can acquire the status of a first 
principle determining the key features of quantum mechanics, i.e. besides the Heisenberg uncertainty 
relations also fully determining the non-local Moyal dynamics (i.e. the interfering trajectories) and 
the standard form of quantum collapse (i.e. the interference between the observer and the system 
during the measurement process). This will be the subject of a future work [31]. 

Finally, we determined the explicit form of the quantum conditions for several different quantum 
phase space ordering prescriptions. For all the studied cases the boundary conditions for the quan- 
tum pure state differential equation are completely specified by the values of Vi2ix) and < p{x) >i2 
(eqs. (49,50)). These two quantities also play an important role in the De Broglie-Bohm formulation 
of quantum mechanics [38]. In fact they also uniquely determine the pure state Dc Broglic-Bohm 
quasi- distribution. This property together with the fact that De Broglie-Bohm theory can be for- 
mulated in phase space in terms of Cohen's quasi-distribution formalism [36, 39, 40] leads to the 
conjecture that the results of this paper can be properly extended to the De Broghe-Bohm rep- 
resentation. Hence, an interesting issue for future research would be to determine the explicit 
form of the pure state quantum condition in the De Broglie-Bohm formulation and show (if possi- 
ble) that its solutions display the well known functional form of the pure state De Broglie-Bohm 
quasi-distributions. 
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